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Abstract 



I discuss several theoretical tools which are useful for analyzing perturbative and 
non-perturbative problems in quantum chromodynamics, including (a) the light-cone 
Fock expansion, (b) the effective charge ay, (c) conformal symmetry, and (d) com- 
mensurate scale relations. Light-cone Fock-state wavefunctions encode the proper- 
ties of a hadron in terms of its fundamental quark and gluon degrees of freedom. 
Given the proton's light-cone wavefunctions, one can compute not only the quark 
and gluon distributions measured in deep inelastic lepton-proton scattering, but also 
the multi-parton correlations which control the distribution of particles in the proton 
fragmentation region and dynamical higher twist effects. Light-cone wavefunctions 
also provide a systematic framework for evaluating exclusive hadronic matrix ele- 
ments, including timelike heavy hadron decay amplitudes and form factors. The ay 
coupling, defined from the QCD heavy quark potential, provides a physical expansion 
parameter for perturbative QCD with an analytic dependence on the fermion masses 
which is now known to two-loop order. Conformal symmetry provides a template for 
QCD predictions, including relations between observables which are present even in 
a theory which is not scale invariant. Commensurate scale relations are perturbative 
QCD predictions based on conformal symmetry relating observable to observable at 
fixed relative scale. Such relations have no renormalization scale or scheme ambiguity. 



1 Introduction 

A primary goal of both high energy and nuclear physics is to unravel the structure and 
dynamics of nucleons and nuclei in terms of their fundamental quark and gluon degrees 
of freedom. Our present empirical knowledge of the quark and gluon distributions of 
the proton has revealed a remarkably complex substructure. It is helpful to categorize 
the parton distributions as "intrinsic" -pertaining to the composition of the target 
hadron, and "extrinsic", reflecting the substructure of the individual quarks and 
gluons themselves. For example, the u[x) and d[x) antiquark distributions of the 
proton at Q'^ ~ 10 GeV^ to be quite different in shape[|I| and thus must reflect 
dynamics intrinsic to the proton's structure. If the sea quarks were generated solely by 
perturbative QCD evolution via gluon splitting, the anti-quark distributions would be 
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isospin symmetric. Evidence for a difference between the s{x) and s{x) distributions 
has also been claimed. ||^ Gluons carry a significant fraction of the proton's spin 
as well as its momentum. Since gluon exchange between valence quarks contributes 
to the p — A mass splitting, it follows that the gluon distributions must be cannot 
be solely accounted for by gluon bremsstrahlung from individual quarks, the process 
responsible for DGLAP evolutions of the structure functions. Similarily, in the case of 
heavy quarks, ss, cc, bb, the diagrams in which the sea quarks are multiply connected 
to the valence quarks are intrinsic to the proton structure itself. Thus neither gluons 
nor sea quarks are solely generated by DGLAP evolution, and one cannot define a 
resolution scale Qq where the sea or gluon degrees of freedom can be neglected. There 
have also been surprises associated with the chirality distributions Aq = qyi — qi/-\ 
of the valence quarks which again show that a simple valence quark approximation 
to nucleon spin structure functions is far from the actual dynamical situation. For a 
recent discussion and references, see Ref. 0. 

A traditional focus of QCD has been on hard inclusive processes and jet physics 
where perturbative methods and leading-twist factorization provide predictions up to 
next-to-next-to leading order (NNLO) with very good precision. More recently, the 
domain of reliable perturbative QCD predictions has been extended to much more 
complex phenomena, such as a fundamental understanding of the hard QCD BFKL 
pomeron in deep inelastic scattering at small x^j and hard diffractive processes, such 
as 7*]9 p^p. In these lectures I will discuss applications of QCD where the non- 
perturbative composition of hadrons in terms of their quark and gluon degrees of 
freedom play a crucial role, for example the x^j-dependence of structure functions 
measured in deep inelastic scattering, exclusive and semi-exclusive processes such 
as form factors, two-photon processes, elastic scattering at fixed 9cm, as well as the 
semi-leptonic decays of heavy hadrons. The analysis of QCD processes at the am- 
plitude level is a challenging problem, mixing issues involving non-perturbative and 
perturbative dynamics. However, a number of tools are available: 

1. The Light-Cone Fock expansion provides a frame-independent representation 
of a hadrons in terms of a set of wavefunctions {ipn/H{xi, k±i, Aj)} describing its com- 
position into relativistic quark and gluon constituents. The light-cone wavefunctions 
can be derived from the eigensolutions of the QCD Hamiltonian defined at fixed light- 
cone time T = t + z/c. Structure functions are obtained from the sum over absolute 
squares of the light-cone wavefunctions. Spacelike form factors and semi-leptonic de- 
cay amplitudes can be written as exact identities in terms of the convolution of the 
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light-cone wavefunctions. 

2. Factorization theorems for hard exclusive, semi-exclusive, and diffractive pro- 
cesses allow a rigorous separation of soft non-perturbative dynamics of the bound 
state hadrons from the hard dynamics of the perturbatively-calculable quark-gluon 
scattering amplitude Tj^\ The key non-perturbative input is the gauge and frame 
independent hadron distribution amplitude [Q defined as the integral over transverse 
momenta of the valence (lowest particle number) Fock wavefunction; e.g. for the pion 

Q)^l d?k^ i^f/M. ku. A) (1) 

where the global cutoff A is identified with the resolution Q. The distribution ampli- 
tude controls leading-twist exclusive amplitudes at high momentum transfer, and it 
can be related to the gauge-invariant Bethe-Salpeter wavefunction at equal light-cone 
time r = . Thus hard exclusive hadronic amplitudes such as quarkonium decay, 
heavy hadron decay, and scattering amplitudes where the hadrons are scattered with 
momentum transfer can be factorized as the convolution of the light-cone Fock state 
wavefunctions with quark-gluon matrix elements 

/n n / n \ / n \ 

n d^ki_ n dxiS 1 - I] a;, 5 iY,k±i] 

n " i=l i=l \ i=l / \i=l / 

xij[%{xiM,,X,)Ti^^ . (2) 

Here is the underlying quark-gluon subprocess scattering amplitude, where the 
(incident or final) hadrons are replaced by quarks and gluons with momenta XiP^, 
XiP± + k±i and invariant mass above the separation scale Ai"^ > A^. 

3. The logarithmic evolution of hadron distribution amplitudes (pn^Xi^Q) can be 
derived from the perturbatively-computable tail of the valence light-cone wavefunc- 
tion in the high transverse momentum regime. [|] 

4. Conformal symmetry provides a template for QCD predictions, leading to 
relations between observables which are present even in a theory which is not scale 
invariant. Thus an important guide in QCD analyses is to identify the underlying 
conformal relations of QCD which are manifest if we drop quark masses and effects 
due to the running of the QCD couplings. In fact, if QCD has an infrared fixed point 
(vanishing of the Cell Mann-Low function at low momenta), the theory will closely 
resemble a scale-free conformally symmetric theory in many applications. 

5. Commensurate scale relations are perturbative QCD predictions which relate 
observable to observable at fixed relative scale, such as the "generalized Crewther 
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relation" , which connects the Bjorken and Gross-Llewellyn Smith deep inelastic scat- 
tering sum rules to measurements of the e^e~ annihilation cross section. The relations 
have no renormalization scale or scheme ambiguity. The coefficients in the perturba- 
tive series for commensurate scale relations are identical to those of conformal QCD; 
thus no infrared renormalons are present. All non-conformal effects are absorbed by 
fixing the ratio of the respective momentum transfer and energy scales. In the case of 
fixed-point theories, commensurate scale relations relate both the ratio of couplings 
and the ratio of scales as the fixed point is approached. 

5. ay Scheme. A natural scheme for defining the QCD coupling in exclusive 
and other processes is the ay scheme defined from heavy quark interactions. All 
vacuum polarization corrections due to fermion pairs are then automatically and 
analytically incorporated into the Gell Mann-Low function, thus avoiding the problem 
of explicitly computing and resumming quark mass corrections related to the running 
of the coupling. 

6. The Abelian Correspondence Principle. One can consider QCD predictions as 
analytic functions of the number of colors Nc and flavors Np- In particular, one can 
show at all orders of perturbation theory that PQCD predictions reduce to those of 
an Abelian theory at Nq — ^ with a = Cpois and Np = Np/TCp held fixed. ||93[] 
There is thus a deep connection between QCD processes and their corresponding 
QED analogs. 

2 The Light-Cone Fock Expansion in QCD 

In a relativistic collision, an incident hadron projectile presents itself as an ensemble 
of coherent states containing various numbers of quark and gluon quanta. Thus when 
a laser beam crosses a proton at fixed "light-cone" time T = t + z/c = x^ + x^, it 
encounters a baryonic state with a given number of quarks, anti-quarks, and gluons 
in flight with fig — rig = 3. The natural formalism for describing these hadronic 
components in QCD is the light-cone Fock representation obtained by quantizing the 
theory at fixed r.0 For example, the proton state has the Fock expansion 
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representing the expansion of the exact QCD eigenstate on a non-interacting quark 
and gluon basis. The probabihty amphtude for each such n-particle state of on-mass 
shell quarks and gluons in a hadron is given by a light-cone Fock state wavefunc- 
tion ipn/H{xi, k±i, Aj), where the constituents have longitudinal light-cone momentum 
fractions 



kt k^ + kf 

Xi 



n 



p+ pO _|_ pz 

relative transverse momentum 



E^^ = l' (4) 



i=l 



k±i , '^k±i = 0± , (5) 



i=l 



and helicities Aj. The effective lifetime of each configuration in the laboratory frame 



is where 



" ^2 _|_ ™2 

mi=j:^^±^<a' (6) 



is the off-shell invariant mass and A is a global ultraviolet regulator. The form of 
'^l!/^Hi^i^k±i, Xi) is invariant under longitudinal boosts; i.e., the light-cone wavefunc- 
tions expressed in the relative coordinates Xj and k±i are independent of the total 
momentum P~^, P± of the hadron. 

The interactions of the proton reflects an average over the interactions of its fluc- 
tuating states. For example, a valence state with small impact separation, and thus 
a small color dipole moment, would be expected to interact weakly in a hadronic or 
nuclear target reflecting its color transparency. The nucleus thus filters differentially 
different hadron components, p, |^ The ensemble {ipn/n} of such light-cone Fock wave- 
functions is a key concept for hadronic physics, providing a conceptual basis for rep- 
resenting physical hadrons (and also nuclei) in terms of their fundamental quark and 
gluon degrees of freedom. Given the V'i/i^) we can construct any spacelike electromag- 
netic or electroweak form factor from the diagonal overlap of the LC wavef unctions. 
Similarly, the matrix elements of the currents that define quark and gluon structure 
functions can be computed from the integrated squares of the LC wavefunctions. 0, |^ 
In general the LC ultraviolet regulators provide a factorization scheme for elastic and 
inelastic scattering, separating the hard dynamical contributions with invariant mass 
squared Ai^ > A"^ from the soft physics with Ai^ < which is incorporated in the 
nonperturbative LC wavefunctions. (Similarly, the DGLAP evolution of quark and 
gluon distributions can be derived by computing the variation of the Fock expansion 
with respect to A^.P]) 
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The light-cone Fock formahsm is derived in the following way: one first constructs 
the light-cone time evolution operator = — and the invariant mass operator 
Hlc = P^P^ — P± in light-cone gauge A'^ = from the QCD Lagrangian. The total 
longitudinal momentum P+ = P^ + P^ and transverse momenta Pi are conserved, i.e. 
are independent of the interactions. The matrix elements of Hlc on the complete or- 
thonormal basis {| n)} of the free theory H^^j = Hicid = 0) can then be constructed. 
The matrix elements {n \ H^c I ^) connect Fock states differing by 0, 1, or 2 quark 
or gluon quanta, and they include the instantaneous quark and gluon contributions 
imposed by eliminating dependent degrees of freedom in light-cone gauge. 

It is thus important to not only compute the spectrum of hadrons and gluonic 
states, but also to determine the wavefunction of each QCD bound state in terms of 
its fundamental quark and gluon degrees of freedom. If we could obtain such nonper- 
turbative solutions of QCD, then we could compute the quark and gluon structure 
functions and distribution amplitudes which control hard-scattering inclusive and 
exclusive reactions as well as calculate the matrix elements of currents which under- 
lie electroweak form factors and the weak decay amplitudes of the light and heavy 
hadrons. The light-cone wavefunctions also determine the multi-parton correlations 
which control the distribution of particles in the proton fragmentation region as well 
as dynamical higher twist effects. Thus one can analyze not only the deep inelastic 
structure functions but also the fragmentation of the spectator system. Knowledge 
of hadron wavefunctions would also open a window to a deeper understanding of the 
physics of QCD at the amplitude level, illuminating exotic effects of the theory such 
as color transparency, intrinsic heavy quark effects, hidden color, diffractive processes, 
and the QCD van der Waals interactions. 

Solving a quantum field theory such as QCD is clearly not easy. However, highly 
nontrivial, one-space one-time relativistic quantum field theories which mimic many 
of the features of QCD, have already been completely solved using light-cone Hamil- 
tonian methods.]^ Virtually any (1+1) quantum field theory can be solved using the 
method of Discretized Light-Cone-Quantization (DLCQ) |jlO|, |ll| where the matrix 
elements (n \ H^]j \ m^, are made discrete in momentum space by imposing periodic 
or anti-periodic boundary conditions in x~ = — x'^ and x^. Upon diagonalization 
of Hici the eigenvalues provide the invariant mass of the bound states and eigen- 
states of the continuum. In DLCQ, the Hamiltonian Hlci which can be constructed 
from the Lagrangian using light-cone time quantization, is completely diagonalized, 
in analogy to Heisenberg's solution of the eigenvalue problem in quantum mechan- 
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ics. The quantum field theory problem is rendered discrete by imposing periodic or 
anti-periodic boundary conditions. The eigenvalues and eigensolutions of collinear 
QCD then give the complete spectrum of hadrons, nuclei, and gluonium and their 
respective light-cone wavefunctions. A beautiful example is "collinear" QCD: a vari- 
ant of QCD{3 + 1) defined by dropping all of interaction terms in H^c^ involving 
transverse momentarily Even though this theory is effectively two-dimensional, the 
transversely-polarized degrees of freedom of the gluon field are retained as two scalar 
fields. Antonuccio and Dalley |jl3l have used DLCQ to solve this theory. The diag- 
onalization of H^c provides not only the complete bound and continuum spectrum 
of the collinear theory, including the gluonium states, but it also yields the complete 
ensemble of light-cone Fock state wavefunctions needed to construct quark and gluon 
structure functions for each bound state. Although the collinear theory is a drastic 
approximation to physical QCD{3 + 1), the phenomenology of its DLCQ solutions 
demonstrate general gauge theory features, such as the peaking of the wavefunctions 
at minimal invariant mass, color coherence and the helicity retention of leading par- 
tons in the polarized structure functions at x — > 1. The solutions of the quantum 
field theory can be obtained for arbitrary coupling strength, flavors, and colors. 

In practice it is essential to introduce an ultraviolet regulator in order to limit the 
total range of {n \ H^c I such as the "global" cutoff in the invariant mass of the 
free Fock state. One can also introduce a "local" cutoff to limit the change in invari- 
ant mass — A^ml < ^loca\ which provides spectator-independent regularization 
of the sub-divergences associated with mass and coupling renormalization. Recently, 
Hiller, McCartor, and I have shown fl^ that the Pauli-Villars method has advantages 
for regulating light-cone quantized Hamitonian theory. We show that Pauli-Villars 
fields satisfying three spectral conditions will regulate the interactions in the ultravio- 
let, while at same time avoiding spectator-dependent renormalization and preserving 
chiral symmetry. Although gauge theories are usually quantized on the light-cone in 
light-cone gauge = 0, it is also possible and interesting to quantize the theory 



in Feynman gauge |T^. Covariant gauges are advantageous since they preserve the 
rotational symmetry of the gauge interactions. 

The natural renormalization scheme for the QCD coupling is av{Q), the effective 
charge defined from the scattering of two infinitely-heavy quark test charges. This is 
discussed in more detail below. The renormalization scale can then be determined 
from the virtuality of the exchanged momentum, as in the BLM and commensu- 



rate scale methods. |jT6|, |T^, |T8|, Similar effective charges have been proposed by 
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Watson pO| and Czarneckiet a/. pT| 



In principle, we could also construct the wavefunctions of QCD(3+1) starting with 
collinear QCD(1+1) solutions by systematic perturbation theory in Aif, where AH 
contains the terms which produce particles at non-zero k^, including the terms linear 
and quadratic in the transverse momenta k j^j which are neglected in the Hamilton 
Hq of collinear QCD. We can write the exact eigensolution of the full Hamiltonian as 

^(3+1) = ^(1+1) + ]^^rrVq~^^^(i+i) ' 

where 

1 1 1 AT/- 1 



can be represented as the continued iteration of the Lippmann Schwinger resolvant. 
Note that the matrix (M^ — Hq)^^ is known to any desired precision from the DLCQ 
solution of collinear QCD. 

3 Electroweak Matrix Elements and Light-Cone 
Wavefunctions 



Dae Sung Hwang and I have recently shown that exclusive semileptonic -B-decay 
amplitudes, such as i? ^ AiV can be evaluated exactly in the light-cone formalism. 
These timelike decay matrix elements require the computation of the diagonal 



matrix element n —>■ n where parton number is conserved, and the off-diagonal n+1 
n — 1 convolution where the current operator annihilates a qq' pair in the initial B 
wavefunction. See Fig. 0. This term is a consequence of the fact that the time-like 
decay = {pe + p^jY > requires a positive light-cone momentum fraction g+ > 0. 
Conversely for space-like currents, one can choose = 0, as in the Drell-Yan-West 
representation of the space-like electromagnetic form factors. However, as can be seen 
from the explicit analysis of the form factor in a perturbation model, the off-diagonal 
convolution can yield a nonzero q^/q'^ limiting form as q^ —>■ 0. This extra term 
appears specifically in the case of "bad" currents such as J~ in which the coupling 
to qq fluctuations in the light-cone wavefunctions are favored. In effect, the g"*" ^ 
limit generates 6{x) contributions as residues of the n + 1 ^ n — 1 contributions. 
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Figure 1: Exact representation of electroweak decays and time- like form factors in 
the light-cone Fock representation. 



The necessity for such "zero mode" 6{x) terms has been noted by Chang, Root and 
Yan,|2l,Burkardt,|2l and Ji and Choi.|25l 



The off-diagonal n + 1 — n — 1 contributions give a new perspective for the 
physics of i?-decays. A semileptonic decay involves not only matrix elements where 
a quark changes flavor, but also a contribution where the leptonic pair is created 
from the annihilation of a qq' pair within the Fock states of the initial B wavefunc- 
tion. The semileptonic decay thus can occur from the annihilation of a nonvalence 
quark-antiquark pair in the initial hadron. This feature will carry over to exclusive 
hadronic 5-decays, such as —>■ tt~D~^. In this case the pion can be produced from 
the coalescence of a du pair emerging from the initial higher particle number Fock 
wavefunction of the B. The D meson is then formed from the remaining quarks after 
the internal exchange of a W boson. 

In principle, a precise evaluation of the hadronic matrix elements needed for B- 
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decays and other exclusive electroweak decay amplitudes requires knowledge of all of 
the light-cone Fock wavefunctions of the initial and final state hadrons. In the case 



of model gauge theories such as QCD(l+l)|g6l or collinear QCD |[T3[ in one-space 
and one-time dimensions, the complete evaluation of the light-cone wavefunction is 
possible for each baryon or meson bound-state using the DLCQ method. It would be 
interesting to use such solutions as a model for physical 5-decays. 

The existence of an exact formalism for electroweak matrix elements gives a basis 
for systematic approximations and a control over neglected terms. For example, one 
can analyze exclusive semileptonic i?-decays which involve hard internal momentum 
transfer using a perturbative QCD formalism patterned after the analysis of form fac- 
tors at large momentum transfer. 0] The hard-scattering analysis proceeds by writing 
each hadronic wavefunction as a sum of soft and hard contributions 

^„ = rn°'\Ml < A') + ^Pl^''^{Ml > A'), (7) 

where is the invariant mass of the partons in the n-particle Fock state and A is 
the separation scale. The high internal momentum contributions to the wavefunction 
i/j^'^'^ can be calculated systematically from QCD perturbation theory by iterating the 
gluon exchange kernel. The contributions from high momentum transfer exchange to 
the 5-decay amplitude can then be written as a convolution of a hard scattering 
quark-gluon scattering amplitude Th with the distribution amplitudes 0(xj,A), the 
valence wavefunctions obtained by integrating the constituent momenta up to the 
separation scale A4n < A < Q. This is the basis for the perturbative hard scatter- 
ing analyses. p7|, p^ , |29| , ^ In the exact analysis, one can identify the hard PQCD 



contribution as well as the soft contribution from the convolution of the light-cone 
wavefunctions. Furthermore, the hard scattering contribution can be systematically 
improved. For example, off-shell effects can be retained in the evaluation of by 
utilizing the exact light-cone energy denominators. 

Given the solution for the hadronic wavefunctions tp^^ with Ai^ < A^, one can 
construct the wavefunction in the hard regime with Ai"^ > A^ using projection opera- 
tor techniques. Q The construction can be done perturbatively in QCD since only high 
invariant mass, far off-shell matrix elements are involved. One can use this method 
to derive the physical properties of the LC wavefunctions and their matrix elements 
at high invariant mass. Since M."^ = J27=i ( ^^t"^ ) ; this method also allows the 



derivation of the asymptotic behavior of light-cone wavefunctions at large k±, which 
in turn leads to predictions for the fall-off of form factors and other exclusive matrix 
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elements at large momentum transfer, such as the quark counting rules for predict- 
ing the nominal power-law fall-off of two-body scattering amplitudes at fixed 6'cm-@] 
The phenomenological successes of these rules can be understood within QCD if the 
coupling av{Q) freezes in a range of relatively small momentum transfer. [lOf 



4 Other Applications of Light-Cone Quantization 
to QCD Phenomenology 

Diffractive vector meson photoproduction. The light-cone Fock wavefunction repre- 
sentation of hadronic amplitudes allows a simple eikonal analysis of diffractive high 
energy processes, such as 'y*{Q^)p — ^ pp, in terms of the virtual photon and the vector 
meson Fock state light-cone wavefunctions convoluted with the gp gp near-forward 



matrix element. |^T[| One can easily show that only small transverse size 6_l ~ 1/Q of 
the vector meson distribution amplitude is involved. The hadronic interactions are 
minimal, and thus the 'j*{Q'^)N — > pN reaction can occur coherently throughout a 
nuclear target in reactions without absorption or shadowing. The 7*y4 V A process 
thus is a laboratory for testing QCD color transparency. [0^] This is discussed further 



in the next section. 

Regge behavior of structure functions. The light-cone wavefunctions ipn/H of a 
hadron are not independent of each other, but rather are coupled via the equations of 



motion. Antonuccio, Dalley and I [^] have used the constraint of finite "mechanical" 
kinetic energy to derive "ladder relations" which interrelate the light-cone wavefunc- 
tions of states differing by one or two gluons. We then use these relations to derive the 
Regge behavior of both the polarized and unpolarized structure functions at x — 0, 
extending Mueller's derivation of the BFKL hard QCD pomeron from the properties 
of heavy quarkonium light-cone wavefunctions at large Nc QCD.|]34 



Structure functions at large x^j. The behavior of structure functions where one 
quark has the entire momentum requires the knowledge of LC wavefunctions with 
X ^ 1 for the struck quark and x — > for the spectators. This is a highly off- 
shell configuration, and thus one can rigorously derive quark-counting and helicity- 
retention rules for the power-law behavior of the polarized and unpolarized quark and 
gluon distributions in the x — > 1 endpoint domain. It is interesting to note that the 
evolution of structure functions is minimal in this domain because the struck quark is 
highly virtual as x — > 1; i.e. the starting point Qq for evolution cannot be held fixed. 
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but must be larger than a scale of order (m^ + k\)/{l — a;).^, ||, |3^ 

Intrinsic gluon and heavy quarks. The main features of the heavy sea quark-pair 
contributions of the Fock state expansion of light hadrons can also be derived from 
perturbative QCD, since M."^ grows with toq. One identifies two contributions to the 
heavy quark sea, the "extrinsic" contributions which correspond to ordinary gluon 
splitting, and the "intrinsic" sea which is multi-connected via gluons to the valence 



quarks. The intrinsic sea is thus sensitive to the hadronic bound state structure. p6 
The maximal contribution of the intrinsic heavy quark occurs at xq ^ Tn±Q / Yl,i ^± 



where m± = \Jm? + k\_\ i.e. at large xq, since this minimizes the invariant mass 
M."^- The measurements of the charm structure function by the EMC experiment are 
consistent with intrinsic charm at large x in the nucleon with a probability of order 
0.6 ± 0.3%. 1^ Similarly, one can distinguish intrinsic gluons which are associated 
with multi-quark interactions and extrinsic gluon contributions associated with quark 
substructure, [p^j One can also use this framework to isolate the physics of the anomaly 
contribution to the Ellis- Jaffe sum rule. 

Materialization of far-off-shell configurations. In a high energy hadronic colli- 
sions, the highly-virtual states of a hadron can be materialized into physical hadrons 
simply by the soft interaction of any of the constituents. [^] Thus a proton state with 
intrinsic charm \ uudcc) can be materialized, producing a. J/ip &i large xp, by the 
interaction of a light-quark in the target. The production occurs on the front-surface 
of a target nucleus, implying an A^/^ J production cross section at large xp , which 
is consistent with experiment, such as Fermilab experiments E772 and E866. 

Rearrangement mechanism in heavy quarkonium decay. It is usually assumed 
that a heavy quarkonium state such as the J/ip always decays to light hadrons via 
the annihilation of its heavy quark constituents to gluons. However, as Karliner 



and I |^0| have recently shown, the transition J/ip pre can also occur by the 
rearrangement of the cc from the J/ ip into the | qqcc) intrinsic charm Fock state of 
the p or vr. On the other hand, the overlap rearrangement integral in the decay 
ip' piT will be suppressed since the intrinsic charm Fock state radial wavefunction 
of the light hadrons will evidently not have nodes in its radial wavefunction. This 
observation gives a natural explanation of the long-standing puzzle why the J/ip 
decays prominently to two-body pseudoscalar- vector final states, whereas the ip' does 
not. 

Asymmetry of intrinsic heavy quark sea. The higher Fock state of the proton 
I uudss) should resemble a | KA) intermediate state, since this minimizes its invari- 
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ant mass M.. In such a state, the strange quark has a higher mean momentum 
fraction x than the s. ||1T|, ^ Similarly, the helicity intrinsic strange quark in 



this configuration will be anti-aligned with the helicity of the nucleon.^, ^ This 
Q ^ Q asymmetry is a striking feature of the intrinsic heavy-quark sea. 

Comover phenomena. Light-cone wavefunctions describe not only the partons 
that interact in a hard subprocess but also the associated partons freed from the 
projectile. The projectile partons which are comoving {i.e., which have similar ra- 
pidity) with final state quarks and gluons can interact strongly producing (a) leading 
particle effects, such as those seen in open charm hadroproduction; (b) suppression 
of quarkonium|[44|| in favor of open heavy hadron production, as seen in the E772 
experiment; (c) changes in color configurations and selection rules in quarkonium 
hadroproduction, as has been emphasized by Hoyer and Peigne.pSj All of these ef- 
fects violate the usual ideas of factorization for inclusive reactions. Further, more than 
one parton from the projectile can enter the hard subprocess, producing dynamical 
higher twist contributions, as seen for example in Drell-Yan experiments. P6|, ^ 

Jet hadronization in light-cone QCD. One of the goals of nonperturbative analysis 
in QCD is to compute jet hadronization from first principles. The DLCQ solutions 
provide a possible method to accomplish this. By inverting the DLCQ solutions, we 
can write the "bare" quark state of the free theory as | go) = I ^) I Qo) where now 
{\n)} are the exact DLCQ eigenstates of H^c, and | go) are the DLCQ projec- 
tions of the eigensolutions. The expansion in automatically infrared and ultraviolet 
regulated if we impose global cutoffs on the DLCQ basis: < AAi"^ < where 
AAI^ = Aifi — (SA^j)^. It would be interesting to study jet hadronization at the 
amplitude level for the existing DLCQ solutions to QCD (1+1) and collinear QCD. 

Hidden Color. The deuteron form factor at high is sensitive to wavefunction 
configurations where all six quarks overlap within an impact separation bj_i < 0{1/Q)] 
the leading power-law fall off predicted by QCD is F^iQ'^) = /(as(<5^))/(<5^)^, where, 
asymptotically, f{as{Q^)) oc q;s((5^)^"''^'''.|^ The derivation of the evolution equation 
for the deuteron distribution amplitude and its leading anomalous dimension 7 is 
given in Ref. In general, the six-quark wavefunction of a deuteron is a mixture of 
five different color-singlet states. The dominant color configuration at large distances 
corresponds to the usual proton-neutron bound state. However at small impact space 
separation, all five Fock color-singlet components eventually acquire equal weight, 
i.e., the deuteron wavefunction evolves to 80% "hidden color." The relatively large 
normalization of the deuteron form factor observed at large points to sizable 
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hidden color contributions. 150 



Spin-Spin Correlations in Nucleon-Nucleon Scattering and the Charm Threshold. 
One of the most striking anomahes in elastic proton-proton scattering is the large spin 



correlation A^vat observed at large angles. [ 51 1 At ~ 5 GeV, the rate for scattering 
with incident proton spins parallel and normal to the scattering plane is four times 
larger than that for scattering with anti-parallel polarization. This strong polarization 
correlation can be attributed to the onset of charm production in the intermediate 



state at this energy. ||52[| The intermediate state \uuduudcc) has odd intrinsic parity 
and couples to the J = S = 1 initial state, thus strongly enhancing scattering when 
the incident projectile and target protons have their spins parallel and normal to the 
scattering plane. The charm threshold can also explain the anomalous change in color 
transparency observed at the same energy in quasi-elastic pp scattering. A crucial 
test is the observation of open charm production near threshold with a cross section 
of order of 1/ib. 



5 Hard Exclusive Reactions 

Exclusive hard-scattering reactions and hard diffractive reactions are now giving a 
valuable window into the structure and dynamics of hadronic amplitudes. Recent 
measurements of the photon-to-pion transition form factor at CLEO,||S^ the diffrac- 
tive dissociation of pions into jets at Fermilab,[^ diffractive vector meson leptopro- 



duction at Fermilab and HERA, and the new program of experiments on exclusive 
proton and deuteron processes at Jefferson Laboratory are now yielding fundamen- 
tal information on hadronic wavef unctions, particularly the distribution amplitude 
of mesons. Such information is also critical for interpreting exclusive heavy hadron 
decays and the matrix elements and amplitudes entering CP- violating processes at 
the B factories. 

There has been much progress analyzing exclusive and diffractive reactions at 
large momentum transfer from first principles in QCD. Rigorous statements can be 
made on the basis of asymptotic freedom and factorization theorems which separate 
the underlying hard quark and gluon subprocess amplitude from the nonperturbative 
physics incorporated into the process-independent hadron distribution amplitudes 
'/'//(^^j, the valence light-cone wavefunctions integrated over < Q"^. An im- 

portant new application is the recent analysis of hard exclusive B decays by Beneke, 



et al.^^ Key features of such analyses are: (a) evolution equations for distribu- 



15 



tion amplitudes which incorporate the operator product expansion, renormahzation 
group invariance, and conformal symmetry; |Q, ^7\, 0, ^ (b) hadron hehc- 



ity conservation which follows from the underlying chiral structure of QCD;[^ (c) 
color transparency, which eliminates corrections to hard exclusive amplitudes from 
initial and final state interactions at leading power and reflects the underlying gauge 
theoretic basis for the strong interactions; P2|, and (d) hidden color degrees of 
freedom in nuclear wavef unctions, which reflects the color structure of hadron and 



nuclear wavef unctions. ||49|] There have also been recent advances eliminating renor- 
mahzation scale ambiguities in hard-scattering amplitudes via commensurate scale 
relations |6^, ^ |6^ which connect the couplings entering exclusive amplitudes to the 



ay coupling which controls the QCD heavy quark potential.]^ The postulate that 
the QCD coupling has an infrared fixed-point can explain the applicability of confor- 
mal scaling and dimensional counting rules to physical QCD processes. [|^, The 
field of analyzable exclusive processes has recently been expanded to a new range of 
QCD processes, such as electroweak decay amplitudes, highly virtual diffractive pro- 
cesses such as 7*p PP,^^ iH and semi-exclusive processes such as 7*p n^X 



[frOl , [7l| , ^ where the 7r+ is produced in isolation at large pt- 

The natural renormahzation scheme for the QCD coupling in hard exclusive pro- 
cesses is av{Q), the effective charge defined from the scattering of two infinitely-heavy 
quark test charges. The renormahzation scale can then be determined from the vir- 
tuality of the exchanged momentum of the gluons, as in the BLM and commensurate 



scale methods. |Tq, 63, 64, 65 



The main features of exclusive processes to leading power in the transferred mo- 
menta are: 

(1) The leading power fall-off is given by dimensional counting rules for the hard- 
scattering amplitude: Th ~ where n is the total number of fields (quarks. 



leptons, or gauge fields) participating in the hard scattering. |]67| , |68[| Thus the reaction 
is dominated by subprocesses and Fock states involving the minimum number of 
interacting fields. The hadronic amplitude follows this fall-off modulo logarithmic 
corrections from the running of the QCD coupling, and the evolution of the hadron 
distribution amplitudes. In some cases, such as large angle pp — > pp scattering, 
pinch contributions from multiple hard-scattering processes must also be included. 
The general success of dimensional counting rules implies that the effective coupling 
av{Q*) controlling the gluon exchange propagators in are frozen in the infrared, 
i. e., have an infrared fixed point, since the effective momentum transfers Q* exchanged 
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by the gluons are often a small fraction of the overall momentum transfer. The 
pinch contributions are suppressed by a factor decreasing faster than a fixed power. p7[| 

(2) The leading power dependence is given by hard-scattering amplitudes Th 
which conserve quark helicity. [0, |7^ Since the convolution of Th with the light-cone 
wavefunctions projects out states with Lz = 0, the leading hadron amplitudes con- 
serve hadron helicity; i.e., the sum of initial and final hadron helicities are conserved. 

(3) Since the convolution of the hard scattering amplitude Th with the light-cone 
wavefunctions projects out the valence states with small impact parameter, the es- 
sential part of the hadron wavefunction entering a hard exclusive amplitude has a 
small color dipole moment. This leads to the absence of initial or final state interac- 
tions among the scattering hadrons as well as the color transparency, of quasi-elastic 
interactions in a nuclear target. [0^, |62| For example, the amplitude for diffractive 
vector meson photoproduction 7*(Q^)p pp, can be written as convolution of the 
virtual photon and the vector meson Fock state light-cone wavefunctions the gp gp 
near-forward matrix element. [^| One can easily show that only small transverse size 
6x ~ l/Q of the vector meson distribution amplitude is involved. The sum over the 
interactions of the exchanged gluons tend to cancel reflecting its small color dipole 
moment. Since the hadronic interactions are minimal, the 'y*{Q'^)N pN reaction 
at large can occur coherently throughout a nuclear target in reactions without ab- 
sorption or final state interactions. The 'y*A V A process thus provides a natural 
framework for testing QCD color transparency. Evidence for color transparency in 



such reactions has been found by Fermilab experiment E665.[75 



Diffractive multi-jet production in heavy nuclei provides a novel way to measure 
the shape of the LC Fock state wavefunctions and test color transparency. For exam- 
ple, consider the reaction ||, 0, ^ ixA Jeti + Jet2 + A! at high energy where the 
nucleus A! is left intact in its ground state. The transverse momenta of the jets have 
to balance so that k±i + k±2 = Q± < R~^A ? and the light-cone longitudinal momentum 
fractions have to add to Xi + X2 ~ 1 so that Ap^ < Ra^- The process can then occur 
coherently in the nucleus. Because of color transparency, i.e., the cancelation of color 
interactions in a small-size color-singlet hadron, the valence wavefunction of the pion 
with small impact separation will penetrate the nucleus with minimal interactions, 
diffracting into jet pairs. P The Xi = x, X2 = I — x dependence of the di-jet distri- 
butions will thus reflect the shape of the pion distribution amplitude; the k_i_i — k±2 
relative transverse momenta of the jets also gives key information on the underlying 
shape of the valence pion wavefunction. |0, ^ The QCD analysis can be confirmed 
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by the observation that the diffractive nuclear amphtude extrapolated to t = is 
linear in nuclear number A, as predicted by QCD color transparency. The integrated 
diffractive rate should scale as A^/R^ ~ A^^^. A diffractive dissociation experiment 
of this type, E791, is now in progress at Fermilab using 500 GeV incident pions on 



nuclear targets. |54| The preliminary results from E791 appear to be consistent with 
color transparency. The momentum fraction distribution of the jets is consistent with 
a valence light-cone wavefunction of the pion consistent with the shape of the asymp- 
totic distribution amplitude, (l)^^^'^^{x) = ^/3fT,x{l — x). As discussed below, data 
from CLEO|Q for the 77* transition form factor also favor a form for the pion 

distribution amplitude close to the asymptotic solution H to the perturbative QCD 



evolution equation. JT^, El) eHj l1' 13 ^^^^ ^^^o be interesting to study diffractive 
tri-jet production using proton beams pA — > Jeti + Jet2 -l- Jets + ^ to determine the 
fundamental shape of the 3-quark structure of the valence light-cone wavefunction 
of the nucleon at small transverse separation. [0 One interesting possibility is that 
the distribution amplitude of the A(1232) for = 1/2,3/2 is close to the asymp- 
totic form X1X2X3, but that the proton distribution amplitude is more complex. This 
would explain why the p — > A transition form factor appears to fall faster at large 
than the elastic p ^ p and the other p — > A^* transition form factors. [|81[ Con- 
versely, one can use incident real and virtual photons: 7* A Jeti + Jet2 + A' to 
confirm the shape of the calculable light-cone wavefunction for transversely-polarized 
and longitudinally-polarized virtual photons. Such experiments will open up a direct 
window on the amplitude structure of hadrons at short distances. 

There are a large number of measured exclusive reactions in which the empiri- 
cal power law fall-off predicted by dimensional counting and PQCD appears to be 
accurate over a large range of momentum transfer. These include processes such as 
the proton form factor, time-like meson pair production in e~^e~ and 77 annihilation, 
large-angle scattering processes such as pion photoproduction 7p t^^P-, and nuclear 
processes such as the deuteron form factor at large momentum transfer and deuteron 
photodisintegration. A spectacular example is the recent measurements at CESR 
of the photon to pion transition form factor in the reaction 67 evr^. |Q As predicted 
by leading twist QCDg] Q'^F^^o{Q'^) is essentially constant for 1 GeV^ < < 10 
GeV^. Further, the normalization is consistent with QCD at NLO if one assumes 
that the pion distribution amplitude takes on the form (p^^"^"^^ (x) = ^/S/t^x^I — x) 
which is the asymptotic solution 0] to the evolution equation for the pion distribution 



amphtude. 77, Bfl, 80 
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The measured deuteron form factor and the deuteron photodisintegration cross 
section appear to follow the leading-twist QCD predictions at large momentum trans- 



fers in the few GeV region. |8^ The normalization of the measured deuteron form 
factor is large compared to model calculations |50] assuming that the deuteron's six- 
quark wavefunction can be represented at short distances with the color structure of 
two color singlet baryons. This provides indirect evidence for the presence of hidden 
color components as required by PQCD.pQ 



If the pion distribution amplitude is close to its asymptotic form, then one can 
predict the normalization of exclusive amplitudes such as the spacelike pion form fac- 
tor Q'^ Ft^{Q'^) . Next-to-leading order predictions are now becoming available which 
incorporate higher order corrections to the pion distribution amplitude as well as the 
hard scattering amplitude. [0, 8^ However, the normalization of the PQCD pre- 



diction for the pion form factor depends directly on the value of the effective coupling 
av{Q*) at momenta Q*^ ^ g2/20. Assuming av{Q*) ^ 0.4, the QCD LO predic- 
tion appears to be smaller by approximately a factor of 2 compared to the presently 
available data extracted from the original pion electroproduction experiments from 



CEA.[^ A definitive comparison will require a careful extrapolation to the pion pole 
and extraction of the longitudinally polarized photon contribution of the ep n^n 
data. 

An important debate has centered on whether processes such as the pion and 
proton form factors and elastic Compton scattering jp — > 7p might be dominated by 



higher twist mechanisms until very large momentum transfers. [p7| , paj , |89| For example, 
if one assumes that the light-cone wavefunction of the pion has the form ^/'softla^, k±) = 
A exp{—b ^^^^^ ) , then the Feynman endpoint contribution to the overlap integral at 
small k± and x ~ 1 will dominate the form factor compared to the hard-scattering 
contribution until very large Q^. However, the above form of ipsoit{x, k±) has no 
suppression at kj_ = for any x; i.e., the wavefunction in the hadron rest frame does 
not fall-off at all for kj_ = and k^ — oo. Thus such wavefunctions do not represent 
well soft QCD contributions. Furthermore, endpoint contributions will be suppressed 
by the QCD Sudakov form factor, reflecting the fact that a near-on-shell quark must 
radiate if it absorbs large momentum. If the endpoint contribution dominates proton 
Compton scattering, then both photons will interact on the same quark line in a 
local fashion and the amplitude is real, in strong contrast to the QCD predictions 
which have a complex phase structure. The perturbative QCD predict ions ||9^ for the 
Compton amplitude phase can be tested in virtual Compton scattering by interference 
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with Bethe-Heitler processes. |^ It should be noted that there is no apparent endpoint 
contribution which could explain the success of dimensional counting in large angle 
pion photoproduction. 

It is interesting to compare the corresponding calculations of form factors of 
bound states in QED. The soft wavefunction is the Schrodinger-Coulomb solution 
tpis{k) oc (1 + / {amj-edY)"'^ , and the full wavefunction, which incorporates trans- 
versely polarized photon exchange, only differs by a factor (1 + p^/rn^^^). Thus 
the leading twist dominance of form factors in QED occurs at relativistic scales 



Q > m^^^.\92\ Furthermore, there are no extra relative factors of a in the hard- 



scattering contribution. If the QCD coupling ay has an infrared fixed point, then the 
fall-off of the valence wavefunctions of hadrons will have analogous power-law forms, 
consistent with the Abelian correspondence principle. |^3|] If power-law wavefunctions 
are indeed applicable to the soft domain of QCD then, the transition to leading-twist 
power law behavior will occur in the nominal hard perturbative QCD domain where 
Q2»(A;i),m2. 



6 Semi- Exclusive Processes: New Probes of Hadron 
Structure 

A new class of hard "semi-exclusive" processes of the form A + B ^ C + Y , have 
been proposed as new probes of QCD.[^, ^ ^ These processes are characterized 
by a large momentum transfer t = {pA — PcY and a large rapidity gap between 
the final state particle C and the inclusive system Y . Here A, B and C can be 
hadrons or (real or virtual) photons. The cross sections for such processes factorize 
in terms of the distribution amplitudes of A and C and the parton distributions in 
the target B. Because of this factorization semi-exclusive reactions provide a novel 
array of generalized currents, which not only give insight into the dynamics of hard 
scattering QCD processes, but also allow experimental access to new combinations of 
the universal quark and gluon distributions. 

QCD scattering amplitude for deeply virtual exclusive processes like Compton 
scattering 7*p — > 7p and meson production 7*p Mp factorizes into a hard subpro- 
cess and soft universal hadronic matrix elements. P^ , pD| pT| For example, consider 
exclusive meson electroproduction such as ep en^n (Fig. ^a). Here one takes (as 
in DIS) the Bjorken limit of large photon virtuality, with xb = Q'^/{2mpv) fixed. 
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Figure 2: (a): Factorization of 7*p vr+n into a skewed parton distribution (SPD), 
a hard scattering H and the pion distribution amphtude 07^- (b): Semi-exclusive 
process 7^*)p ^ tt"*"!^. The d-quark produced in the hard scattering H hadronizes 
independently of the spectator partons in the proton, (c): Diagram for the cross sec- 
tion of a generic semi-exclusive process. It involves a hard scattering H, distribution 
amplitudes (pA and (pc and a parton distribution (PD) in the target B. 

while the momentum transfer t = [pp — remains small. These processes involve 
'skewed' parton distributions, which are generalizations of the usual parton distribu- 
tions measured in DIS. The skewed distribution in Fig. |^ describes the emission of 
a M-quark from the proton target together with the formation of the final neutron 
from the d-quark and the proton remnants. As the subenergy s of the scattering 
process 7*^ n^d is not fixed, the amplitude involves an integral over the w-quark 
momentum fraction x. 

An essential condition for the factorization of the deeply virtual meson production 
amplitude of Fig. ^ is the existence of a large rapidity gap between the produced 
meson and the neutron. This factorization remains valid if the neutron is replaced 
with a hadronic system Y of invariant mass My -C W^, where W is the cm. energy 
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of the 7*p process. For My ^ rrip the momentum k' of the c?-quark in Fig. ^ is large 
with respect to the proton remnants, and hence it forms a jet. This jet hadronizes 
independently of the other particles in the final state if it is not in the direction of the 
meson, i.e., if the meson has a large transverse momentum q'j_ = with respect to 
the photon direction in the 7*p cm. Then the cross section for an inclusive system 
Y can be calculated as in DIS, by treating the (i-quark as a final state particle. 

The large furthermore allows only transversally compact configurations of the 
projectile A to couple to the hard subprocess H of Fig. ^b, as it does in exclusive 
processes. Hence the above discussion applies not only to incoming virtual photons 
at large Q^, but also to real photons {Q^ = 0) and in fact to any hadron projectile. 

Let us then consider the general process A + B ^ C + Y , where B and C are 
hadrons or real photons, while the projectile A can also be a virtual photon. In the 
semi-exclusive kinematic limit A.qq.j^, Mg, Mq <^ My, A^ ^ W"^ we have a large 
rapidity gap \yc — yd\ = log ^^a^^a between C and the parton d produced in the hard 
scattering (see Fig. ^c). The cross section then factorizes into the form 

.(^A + B^C + Y) 



dt dxs 

= Y.h/B{xs,f-i')^{Ab-.Cd), (8) 

where t = {q — q'Y and fb/B{xs, /i^) denotes the distribution of quarks, antiquarks and 
gluons b in the target B. The momentum fraction of the struck parton b is fixed 
by kinematics to the value xs = and the factorization scale //^ is characteristic 

of the hard subprocess Ab Cd. 

It is conceptually helpful to regard the hard scattering amplitude H in Fig. |^c as 
a generalized current of momentum q — q' = Pa~ Pci which interacts with the target 
parton b. For A = 7* we obtain a close analogy to standard DIS when particle C 
is removed. With g' ^ we thus find — t — > Q^, My W^, and see that xs goes 
over into xb = Q"^ /{W"^ + Q"^). The possibility to control the value of q' (and hence 
the momentum fraction X5 of the struck parton) as well as the quantum numbers of 
particles A and C should make semi-exclusive processes a versatile tool for studying 
hadron structure. The cross section further depends on the distribution amplitudes 
0c (c-/ Fig. ^), allowing new ways of measuring these quantities. 
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7 Conformal Symmetry as a Template 



Testing quantum chromodynamics to high precision is not easy. Even in processes 
involving high momentum transfer, perturbative QCD predictions are comphcated 
by questions of the convergence of the series, particularly by the presence of "renor- 
malon" terms which grow as n!, reflecting the uncertainty in the analytic form of 
the QCD coupling at low scales. Virtually all QCD processes are complicated by 
the presence of dynamical higher twist effects, including power-law suppressed contri- 
butions due to multi-parton correlations, intrinsic transverse momentum, and finite 
quark masses. Many of these effects are inherently nonperturbative in nature and 
require knowledge of hadron wavefunction themselves. The problem of interpreting 
perturbative QCD predictions is further compounded by theoretical ambiguities due 
to the apparent freedom in the choice of renormalization schemes, renormalizations 
scales, and factorization procedures. 

A central principle of renormalization theory is that predictions which relate phys- 
ical observables to each other cannot depend on theoretical conventions. For example, 
one can use any renormalization scheme, such as the modified minimal subtraction 
dimensional regularization scheme, and any choice of renormalization scale /i to com- 
pute perturbative series observables A and B. However, all traces of the choices of 
the renormalization scheme and scale must disappear when we algebraically eliminate 
the 

'^M5(/^) ^"^^ directly relate A to B. This is the principle underlying "commen- 
surate scale relations" (CSR) [|l7l, which are general leading-twist QCD predictions 
relating physical observables to each other. For example, the "generalized Crewther 
relation" , which is discussed in more detail below, provides a scheme-independent re- 
lation between the QCD corrections to the Bjorken (or Gross Llewellyn-Smith) sum 
rule for deep inelastic lepton-nucleon scattering, at a given momentum transfer Q, to 
the radiative corrections to the annihilation cross section cr e+e-^hadmnsis), at a cor- 
responding "commensurate" energy scale ^/s. |T7| , [95[1 The specific relation between 
the physical scales Q and ^/s reflects the fact that the radiative corrections to each 
process have distinct quark mass thresholds. 

Any perturbatively calculable physical quantity can be used to define an effective 
charge [P^, P7| , PB[ by incorporating the entire radiative correction into its definition. 
For example, the e~^e~'j* —>■ hadrons annihilation to muon pair cross section ratio can 
be written 

i?e+e-(s) = i2°+e-(^)[l + ^], (9) 

IT 
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where R^+e- is the prediction at Born level. Similarly, we can define the entire 
radiative correction to the Bjorken sum rule as the effective charge agj^{Q'^) where Q 
is the corresponding momentum transfer: 



/ dx 
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(10) 

By convention, each effective charge is normalized to as in the weak coupling limit. 
One can define effective charges for virtually any quantity calculable in perturbative 
QCD; e.^'.moments of structure functions, ratios of form factors, jet observables, and 
the effective potential between massive quarks. In the case of decay constants of the 
Z or the r, the mass of the decaying system serves as the physical scale in the effective 
charge. In the case of multi-scale observables, such as the two-jet fraction in e'^e~ 
annihilation, the multiple arguments of the effective coupling a2jetis, y) correspond 
to the overall available energy s variables such as y — maxij(pi -|- pj)^/s representing 
the maximum jet mass fraction. 

Commensurate scale relations take the general form 

aA{QA)=CAB[aB{QB)] . (11) 

The function Cab(q;b) relates the observables A and B in the conformal limit; i.e., 
Cab gives the functional dependence between the effective charges which would be 
obtained if the theory had zero /3 function. The conformal coefficients can be distin- 
guished from the terms associated with the /3 function at each order in perturbation 
theory from their color and fiavor dependence, or by an expansion about a fixed point. 

The ratio of commensurate scales is determined by the requirement that all terms 
involving the /3 function are incorporated into the arguments of the running cou- 
phngs, as in the original BLM procedure. Physically, the ratio of scales corresponds 
to the fact that the physical observables have different quark threshold and distinct 
sensitivities to fermion loops. More generally, the differing scales are in effect rela- 
tions between mean values of the physical scales which appear in loop integrations. 
Commensurate scale relations are transitive; i.e., given the relation between effective 
charges for observables A and C and C and B, the resulting between A and B is 
independent of C. In particular, transitivity implies A^^ = Aac x Aqb- 

One can consider QCD predictions as functions of analytic variables of the number 
of colors Nc and flavors Np. For example, one can show at all orders of perturbation 
theory that PQCD predictions reduce to those of an Abelian theory at Nc — > with 
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S = CpOLs and iV^ = Nf/TCp held fixed. In particular, CSRs obey the "Abelian 
correspondence principle" in that they give the correct Abehan relations at Nc 0. 

Similarly, commensurate scale relations obey the "conformal correspondence prin- 
ciple" : the CSRs reduce to correct conformal relations when Nc and Np are tuned 
to produce zero (5 function. Thus conformal symmetry provides a template for QCD 
predictions, providing relations between observables which are present even in theo- 
ries which are not scale invariant. All effects of the nonzero beta function are encoded 
in the appropriate choice of relative scales A^^ = Qa/Qb- 

The scale Q which enters a given effective charge corresponds to a physical momen- 
tum scale. The total logarithmic derivative of each effective charge effective charge 
OiA{Q) with respect to its physical scale is given by the Cell Mann-Low equation: 

^.r ' ' = *A(«A(g,m),g/m), (12) 

where the functional dependence of a is specific to its own effective charge. Here m 
refers to the quark's pole mass. The pole mass is universal in that it does not depend 
on the choice of effective charge. The Cell Mann-Low relation is reflexive in that -0^ 
depends on only on the coupling a a at the same scale. It should be emphasized that 
the Cell Mann-Low equation deals with physical quantities and is independent of 
the renormalization procedure and choice of renormalization scale. A central feature 
of quantum chromodynamics is asymptotic freedom; i.e., the monotonic decrease of 
the QCD couphng Q;^(/i^) at large spacelike scales. The empirical test of asymptotic 
freedom is the verification of the negative sign of the Cell Mann-Low function at large 
momentum transfer, which must be true for any effective charge. 
In perturbation theory, 

*, = -v,?)!^i-ri"4-'/'?4+--- (13) 

At large scales ^ rn? , the first two terms are universal and identical to the first 
two terms of the (3 function ip^^^ = = ^^y^ ~ ^^f,'^a^ = Pi, whereas iIj^^ for 
n > 2 is process dependent. The quark mass dependence of the iJj function is analytic, 
and in the case of ay scheme is known to two loops. 

The commensurate scale relation between aA and implies an elegant rela- 
tion between their conformal dependence Cab and their respective Cell Mann Low 
functions: 

^B = ^X ^Pa. (14) 
aaA 
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Thus given the resuh for Npyim/Q) in ay scheme we can use the CSR to derive 
NF,A{f^/Q) for any other effective charge, at least to two loops. The above relation 
also shows that if one effective charge has a fixed point ipA[(y-A{QA^)] = 0, then 
all effective charges B have a corresponding fixed point ^/'^[^^(Q^'^)] = at the 
corresponding commensurate scale and value of effective charge. 

In quantum electrodynamics, the running coupling dQEoiQ'^)^ defined from the 
Coulomb scattering of two infinitely heavy test charges at the momentum transfer t = 
—Q^, is taken as the standard observable. Is there a preferred effective charge which 
we should use to characterize the coupling strength in QCD? In the case of QCD, 
the heavy-quark potential VlQ"^) is defined via a Wilson loop from the interaction 
energy of infinitely heavy quark and antiquark at momentum transfer t = —Q^. 
The relation V{Q'^) = —4:nCFav{Q'^)/Q'^ then defines the effective charge av{Q)- 
As in the corresponding case of Abelian QED, the scale Q of the coupling av{Q) 
is identified with the exchanged momentum. Thus there is never any ambiguity in 
the interpretation of the scale. All vacuum polarization corrections due to fermion 
pairs are incorporated in ay through the usual vacuum polarization kernels which 
depend on the physical mass thresholds. Other observables could be used to define 
the standard QCD coupling, such as the effective charge defined from heavy quark 
radiation. p9[ 

Commensurate scale relations between ay and the QCD radiative corrections to 
other observables have no scale or scheme ambiguity, even in multiple-scale problems 
such as multi-jet production. As is the case in QED, the momentum scale which 
appears as the argument of ay reflect the mean virtuality of the exchanged gluons. 
Furthermore, we can write a commensurate scale relation between ay and an analytic 
extension of the coupling, thus transferring all of the unambiguous scale-fixing 
and analytic properties of the physical ay scheme to the MS coupling. 

An elegant example is the relation between the rate for semi-leptonic -B-decay and 

ay: 



^ _ ^_^^ ay{0.16M,) _ ^^^^ ay{0.16M,) 



TT TT 



(15) 



where Mb is the scheme independent 6— quark pole mass. The coefficient of a|^^(/i) 
in the usual expansion with fi = ruf, is 26.8. 
Some other examples of CSR's at NLO: 

a^i(0.5v/i) , «gi(0.5v/i) 
aR{y/s) = ag,{0.5y/s) ^ 5 (16) 
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an(.V-s) = avd.Syi) + 2.08^^^^:?^^ - T.ie^^MiMl (17) 

TT TT^ 

a^(y^) = a^(0.8Vs) + 2.08-5^^ — - 7.16 ^ (18) 

a,,(v/5) = av(0.8Q) + 1.08^i<2:M _ lo^3^i<2^« (19) 

For numerical purposes in each case we have used Np = 5 and ay = 0.1 to compute 
the NLO correction to the CSR scale. 

Commensurate scale relations thus provide fundamental and precise scheme-independent 
tests of QCD, predicting how observables track not only in relative normalization, but 
also in their commensurate scale dependence. 



8 The Generalized Crewther Relation 



The generalized Crewther relation can be derived by calculating the QCD radiative 
corrections to the deep inelastic sum rules and Re+e- in a convenient renormalization 
scheme such as the modified minimal subtraction scheme MS. One then algebraically 



eliminates ajfg{^). Finally, BLM scale-sett ing[^] is used to eliminate the /^-function 
dependence of the coefficients. The form of the resulting relation between the ob- 
servables thus matches the result which would have been obtained had QCD been a 
conformal theory with zero f3 function. The final result relating the observables is 
independent of the choice of intermediate MS renormalization scheme. 

More specifically, consider the Adler function ||100|| for the e~^e~ annihilation cross 
section 



D{Q') = -UTi'Q'^^niQ'), n(g 



[s)ds 



127r2 74m2 S(s + g2 



(20) 



The entire radiative correction to this function is defined as the effective charge 



D(Q^/^i\a,{^?)) = d(i,«s(Q' 



(21) 
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where Cp = ^2N^ ■ The coefficient Cl(Q^) appears at the third order in perturbation 
theory and is related to the "hght-by-hght scattering type" diagrams. (Hereafter 
will denote the MS scheme strong coupling constant.) 

It is straightforward to algebraically relate ag^iQ"^) to aoiQ"^) using the known 
expressions to three loops in the MS scheme. If one chooses the renormalization scale 
to resum all of the quark and gluon vacuum polarization corrections into aoiQ'^)-, 
then the final result turns out to be remarkably simple (S = 3/4 Ci? a/vr) : 

agAQ) = Sd(Q*) - aliQ*) + aliQ*) + ■■■, (22) 



where 




§ + f C(3) - 16C^(3)) Po 
" l^"-^^^^ " 1^^^ " ^^fC(3) + 80CfC(5) 




(23) 



where in QCD, Ca = Nc = 3 and Cp = 4/3. This relation shows how the coefficient 
functions for these two different processes are related to each other at their respective 
commensurate scales. We emphasize that the MS renormalization scheme is used only 
for calculational convenience; it serves simply as an intermediary between observables. 
The renormalization group ensures that the forms of the CSR relations in perturbative 
QCD are independent of the choice of an intermediate renormalization scheme. 

The Crewther relation was originally derived assuming that the theory is confor- 
mally invariant; i. e., for zero (3 function. In the physical case, where the QCD coupling 
runs, all non-conformal effects are resummed into the energy and momentum transfer 
scales of the effective couplings and agi. The general relation between these two 
effective charges for non-conformal theory thus takes the form of a geometric series 

l-ag, = [l+Sz)(Q*)]"' . (24) 

We have dropped the small light-by-light scattering contributions. This is again a 
special advantage of relating observable to observable. The coefficients are indepen- 
dent of color and are the same in Abelian, non-Abelian, and conformal gauge theory. 
The non-Abelian structure of the theory is reflected in the expression for the scale 
Q*. 

Is experiment consistent with the generalized Crewther relation? Fits ||101|| to the 
experimental measurements of the i?-ratio above the thresholds for the production 
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of cc bound states provide the empirical constraint: aR{^/s = 5.0 GeV)/7r ~ 0.08 ± 
0.03. The prediction for the effective couphng for the deep inelastic sum rules at 
the commensurate momentum transfer Q is then ^^^(Q = 12.33 ± 1.20 GeV)/7r ~ 
aGLs(<5 = 12.33±1.20 GeV)/7r ~ 0.074±0.026. Measurements of the Gross-Llewellyn 
Smith sum rule have so far only been carried out at relatively small values of Q^ |p.02| , 
[L03|| ; however, one can use the results of the theoretical extrapolation | TU^ ] of the 
experimental data presented inlTol: aGLs'^°^(<5 = 12.25 GeV)/7r ^ 0.093 ± 0.042. 
This range overlaps with the prediction from the generalized Crewther relation. It is 
clearly important to have higher precision measurements to fully test this fundamental 
QCD prediction. 



9 General Form of Commensurate Scale Relations 

In general, commensurate scale relations connecting the effective charges for observ- 
ables A and B have the form 



o:a{Qa) = o:b{Qb) I 1 + r^j^ h 



rf,s^^ +■■■], (25) 



A/B 



TT 



where the coefficients r^^^ are identical to the coefficients obtained in a conformally 
invariant theory with /SsiaB) = {d/dlnQ'^)aB{Q'^) = 0. The ratio of the scales 
Qa/Qb is thus fixed by the requirement that the couplings sum all of the effects of 
the non-zero (3 function. In practice the NLO and NNLO coefficients and relative 
scales can be identified from the flavor dependence of the perturbative series; i.e. by 
shifting scales such that the A?"^- dependence associated with Pq = 11 /3C a — 4:/3TpNp 
and Pi = — 34/3C^ + ^CaTeNp + ACpTpNp does not appear in the coefficients. 
Here Ca = Nc, Cp = {N^ - l)/2Nc and Tp = 1/2. The shift in scales which gives 
conformal coefficients in effect pre-sums the large and strongly divergent terms in the 
PQCD series which grow as n!(/?octs)"; ^-e., the infrared renormalons associated with 
coupling-constant renormalization. ||106| , ^ [lOT] , |108| 



The renormalization scales Q* in the BLM method are physical in the sense that 
they reflect the mean virtuality of the gluon propagators. This scale-flxing proce- 
dure is consistent with scale flxing in QED, in agreement with in the Abelian limit, 
Nc — * 0.|T^, ^ |109| , |11U| , |111| The ratio of scales Xa/b = Qa/Qb guarantees that 



the observables A and B pass through new quark thresholds at the same physical 
scale. One can also show that the commensurate scales satisfy the transitivity rule 
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Aa/b = ^A/c^c/B, which ensures that predictions are independent of the choice of an 
intermediate renormahzation scheme or intermediate observable C. 

In general, we can write the relation between any two effective charges at arbi- 
trary scales and /i^ as a correction to the corresponding relation obtained in a 
conformally invariant theory: 

aA{^^A) = CABlasilJ'B)] + pB[aB{lJ'B)]FAB[aB{lJ'B)] (26) 

where 

CabK] =aB + E41«B (27) 

n=l 

is the functional relation when PbIo^b] = 0. In fact, if as approaches a fixed point 
aB where Pb\(^b] = 0, then a a tends to a fixed point given by 

aA^aA = Cab [as]. (28) 

The commensurate scale relation for observables A and B has a similar form, but in 
this case the relative scales are fixed such that the non-conformal term Fab is zero. 
Thus the commensurate scale relation q;a(<5a) = CABldB^QB)] at general commen- 
surate scales is also the relation connecting the values of the fixed points for any 
two effective charges or schemes. Furthermore, as (3 —>■ 0, the ratio of commensurate 
scales Qa/Qb becomes the ratio of fixed point scales Qa/Qb as one approaches the 
fixed point regime. 



10 Implementation of ay Scheme 

The effective charge ay(Q) provides a physically-based alternative to the usual mod- 
ified minimal subtraction (MS) scheme. All vacuum polarization corrections due to 
fermion pairs are incorporated in ay through the usual vacuum polarization kernels 
which depend on the physical mass thresholds. When continued to time-like mo- 
menta, the coupling has the correct analytic dependence dictated by the production 
thresholds in the crossed channel. Since ay incorporates quark mass effects exactly, 
it avoids the problem of explicitly computing and resumming quark mass corrections 
which are related to the running of the coupling. Thus the effective number of fla- 
vors Np{Q/m) is an analytic function of the scale Q and the quark masses m. The 
effects of finite quark mass corrections on the running of the strong coupling were 



first considered by De Riijula and Georgi 112| within the momentum subtraction 
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schemes (MOM) (see also Georgi and Politzer [ |113|| , Shirkov and collaborators ||114 
and Chyla P^ ). 

One important advantage of the physical charge approach is its inherent gauge 
invariance to all orders in perturbation theory. This feature is not manifest in massive 
/^-functions defined in non-physical schemes such as the MOM schemes. A second, 
more practical, advantage is the automatic decoupling of heavy quarks according to 



the Appelquist-Carazzone theorem |p.l6| . 

By employing the commensurate scale relations other physical observables can be 
expressed in terms of the analytic coupling ay without scale or scheme ambiguity. 
This way the quark mass threshold effects in the running of the coupling are taken 
into account by utilizing the mass dependence of the physical ay scheme. In effect, 
quark thresholds are treated analytically to all orders in m? /Q"^; i.e., the evolution of 
the physical ay coupling in the intermediate regions reflects the actual mass depen- 
dence of a physical effective charge and the analytic properties of particle production. 
Furthermore, the definiteness of the dependence in the quark masses automatically 
constrains the scale Q in the argument of the coupling. There is thus no scale ambi- 
guity in perturbative expansions in ay- 

In the conventional MS scheme, the coupling is independent of the quark masses 
since the quarks are treated as either massless or infinitely heavy with respect to the 
running of the coupling. Thus one formulates different effective theories depending 
on the effective number of quarks which is governed by the scale Q; the massless 
/3-function is used to describe the running in between the flavor thresholds. These 
different theories are then matched to each other by imposing matching conditions at 
the scale where the effective number of flavors is changed (normally the quark masses). 
The dependence on the matching scale can be made arbitrarily small by calculating 
the matching conditions to high enough order. For physical observables one can then 
include the effects of finite quark masses by making a higher- twist expansion in iv? / 
and Q"^ /m? for light and heavy quarks, respectively. These higher-twist contributions 
have to be calculated for each observables separately, so that in principle one requires 
an all-order resummation of the mass corrections to the effective Lagrangian to give 
correct results. 

The specification of the coupling and renormalization scheme also depends on the 
definition of the quark mass. In contrast to QED where the on-shell mass provides a 
natural definition of lepton masses, an on-shell definition for quark masses is compli- 
cated by the confinement property of QCD. In this paper we will use the pole mass 
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m(j9^ = m^) = m which has the advantage of being scheme and renormahzation-scale 
invariant. Furthermore, when combined with the ay scheme, the pole mass gives 
predictions which are free of the leading renormalon ambiguity. 

A technical complication of massive schemes is that one cannot easily obtain 
analytic solutions of renormalization group equations to the massive (3 function, and 
the Gell-Mann Low function is scheme-dependent even at lowest order. 

In a recent paper we have presented a two-loop analytic extension of the ay- 
scheme based on the recent results of Ref. ||117|| . The mass effects are in principle 
treated exactly to two-loop order and are only limited in practice by the uncertainties 
from numerical integration. The desired features of gauge invariance and decoupling 
are manifest in the form of the two-loop Gell-Mann Low function, and we give a simple 
fitting-function which interpolates smoothly the exact two-loop results obtained by 



using the adoptive Monte Carlo integrator VEGAS ||118|| . Strong consistency checks 
of the results are performed by comparing the Abelian limit to the well known QED 
results in the on-shell scheme. In addition, the massless as well as the decoupling 
limit are reproduced exactly, and the two-loop Gell-Mann Low function is shown to 
be renormalization scale independent. 

The results of our numerical calculation of Npy in the l^-scheme for QCD and 
QED are shown in Fig. |^. The decoupling of heavy quarks becomes manifest at small 
Q/m, and the massless limit is attained for large Q/m. The QCD form actually 
becomes negative at moderate values of Q/m, a novel feature of the anti-screening 
non-Abelian contributions. This property is also present in the (gauge dependent) 
MOM results. In contrast, in Abelian QED the two-loop contribution to the effective 
number of flavors becomes larger than 1 at intermediate values of Q/m. We also 
display the one-loop contribution Npy (^^^ which monotonically interpolates between 
the decoupling and massless limits. The solid curves displayed in Fig. ^ shows that the 
parameterizations which we used for fltting the numerical results are quite accurate. 



The relation of ay (Q^) to the conventional MS coupling is now known to NNLO,|]TT^ 
but for clarity in this section only the NLO relation will be used. The commensurate 



scale relation is given by 1 120 



O 71 

= ay(g*) + 2^M^, (29) 

TT 

which is valid for Q^ ^ m^. The coefficients in the perturbation expansion have their 



32 




Figure 3: The numerical results for the gauge-invariant Npy in QED (open circles) 
and QCD (triangles) with the best fits superimposed respectively. The dashed 
line shows the one-loop Npy function . For comparison we also show the gauge 
dependent two-loop result obtained in MOM schemes (dash-dot) ||121| , |122|| . At large 
^ the theory becomes effectively massless, and both schemes agree as expected. The 
figure also illustrates the decoupling of heavy quarks at small ^ . 

conformal values, i.e., the same coefficients would occur even if the theory had been 
conformally invariant with f3 = 0. The commensurate scale is given by 

"51 



Q* = Q exp 



6 



(30) 



The scale in the MS scheme is thus a factor ~ 0.4 smaller than the physical scale. 
The coefficient 2Nc/3 in the NLO coefficient is a feature of the non-Abelian couplings 
of QCD; the same coefficient occurs even if the theory were conformally invariant with 
/5o = 0. 

Using the above QCD results, we can transform any NLO prediction given in 



MS scheme to a scale-fixed expansion in ay(Q). We can also derive the connection 
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between the MS and ay schemes for Abehan perturbation theory using the hmit 
Nc^Q with Cpcis and Np/Cp held fixed. Q 

The use of ay and related physically defined effective charges such as ap (to 
NLO the effective charge defined from the (1,1) plaquette, ap is the same as ay) as 
expansion parameters has been found to be valuable in lattice gauge theory, greatly 
increasing the convergence of perturbative expansions relative to those using the bare 



lattice coupling. ||109|| Recent lattice calculations of the T- spectrum [ |123|| have been 
used with BLM scale-fixing to determine a NLO normalization of the static heavy 
quark potential: (8.2GeV) = 0.196(3) where the effective number of light fiavors 
is n/ = 3. The corresponding modified minimal subtraction coupling evolved to the 
Z mass and five fiavors is a^-^{Mz) = 0.1174(24). Thus a high precision value for 
aviQ"^) at a specific scale is available from lattice gauge theory. Predictions for other 
QCD observables can be directly referenced to this value without the scale or scheme 
ambiguities, thus greatly increasing the precision of QCD tests. 

One can also use ay to characterize the coupling which appears in the hard scat- 
tering contributions of exclusive process amplitudes at large momentum transfer, 
such as elastic hadronic form factors, the photon-to-pion transition form factor at 
large momentum transfer HTBI, |19| and exclusive weak decays of heavy hadrons. [|124|] 
Each gluon propagator with four-momentum k^^ in the hard-scattering quark-gluon 
scattering amplitude can be associated with the coupling avik"^) since the gluon 
exchange propagators closely resembles the interactions encoded in the effective po- 
tential VlQ"^). [In Abelian theory this is exact.] Commensurate scale relations can 
then be established which connect the hard-scattering subprocess amplitudes which 
control exclusive processes to other QCD observables. 

We can anticipate that eventually nonperturbative methods such as lattice gauge 
theory or discretized light-cone quantization will provide a complete form for the 
heavy quark potential in QCD. It is reasonable to assume that av{Q) will not diverge 
at small space-like momenta. One possibility is that ay stays relatively constant 
av{Q) — 0.4 at low momenta, consistent with fixed-point behavior. There is, in 
fact, empirical evidence for freezing of the ay coupling from the observed systematic 
dimensional scaling behavior of exclusive reactions. [T^ If this is in fact the case, then 
the range of QCD predictions can be extended to quite low momentum scales, a 
regime normally avoided because of the apparent singular structure of perturbative 
extrapolations. 

There are a number of other advantages of the \^-scheme: 
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1. Perturbative expansions in ay with the scale set by the momentum transfer 
cannot have any /3-function dependence in their coefficients since all running 
coupling effects are already summed into the definition of the potential. Since 
coefficients involving i3q cannot occur in an expansions in ay, the divergent 
infrared renormalon series of the form ay/S^nl cannot occur. The general con- 
vergence properties of the scale Q* as an expansion in ay is not known. 

2. The effective coupling aviQ"^) incorporates vacuum polarization contributions 
with finite fermion masses. When continued to time-like momenta, the coupling 
has the correct analytic dependence dictated by the production thresholds in 
the t channel. Since ay incorporates quark mass effects exactly, it avoids the 
problem of explicitly computing and resumming quark mass corrections. 



3. The ay coupling is the natural expansion parameter for processes involving non- 
relativistic momenta, such as heavy quark production at threshold where the 
Coulomb interactions, which are enhanced at low relative velocity v as nay/v, 
need to be re-summed. ||125| , |126| , p.27|| The effective Hamiltonian for nonrelativis- 



tic QCD is thus most naturally written in ay scheme. The threshold corrections 
to heavy quark production in e~^e~ annihilation depend on ay at specific scales 
Q*. Two distinct ranges of scales arise as arguments of ay near threshold: the 
relative momentum of the quarks governing the soft gluon exchange responsible 
for the Coulomb potential, and a high momentum scale, induced by hard gluon 



exchange, approximately equal to twice the quark mass for the corrections. [|12(: 
One thus can use threshold production to obtain a direct determination of ay 
even at low scales. The corresponding QED results for r pair production allow 
for a measurement of the magnetic moment of the r and could be tested at a 
future r-charm factory. | |125| , |126 



We also note that computations in different sectors of the Standard Model have 
been traditionally carried out using different renormalization schemes. However, in 
a grand unified theory, the forces between all of the particles in the fundamental 
representation should become universal above the grand unification scale. Thus it 
is natural to use ay as the effective charge for all sectors of a grand unified theory, 
rather than in a convention-dependent coupling such as ajjg. 



35 



11 The Analytic Extension of the MS Scheme 



The standard MS scheme is not an analytic function of the renormahzation scale at 
heavy quark thresholds; in the running of the coupling the quarks are taken as mass- 
less, and at each quark threshold the value of Np which appears in the (3 function is 
incremented. Thus Eq. (^) is technically only valid far above a heavy quark thresh- 
old. However, we can use this commensurate scale relation to define an extended 



MS scheme which is continuous and analytic at any scale. The new modified scheme 
inherits all of the good properties of the ay scheme, including its correct analytic 
properties as a function of the quark masses and its unambiguous scale fixing. [ p^20|] 
Thus we define 

otME\Q) = av{Q > + -T, \ , (31) 

O TT 



for all scales Q. This equation not only provides an analytic extension of the MS 
and similar schemes, but it also ties down the renormahzation scale to the physical 
masses of the quarks as they enter into the vacuum polarization contributions to ay- 



The modified scheme Oms provides an analytic interpolation of conventional MS 
expressions by utilizing the mass dependence of the physical ay scheme. In effect, 
quark thresholds are treated analytically to all orders in m? /Q"^; i.e., the evolution of 
the analytically extended coupling in the intermediate regions reflects the actual mass 
dependence of a physical effective charge and the analytic properties of particle pro- 
duction. Just as in Abelian QED, the mass dependence of the effective potential and 
the analytically extended scheme Oms reflects the analyticity of the physical thresh- 
olds for particle production in the crossed channel. Furthermore, the definiteness 
of the dependence in the quark masses automatically constrains the renormahzation 
scale. There is thus no scale ambiguity in perturbative expansions in ay or a^g. 

In leading order the effective number of flavors in the modified scheme is given 



to a very good approximation by the simple form ||120" | 

. (32) 




Thus the contribution from one flavor is ^ 0.5 when the scale Q equals the quark 
mass rrii. The standard procedure of matching aMs(/i) at the quark masses serves as 
a zeroth-order approximation to the continuous Np. 

Adding all flavors together gives the total A^^°^(Q) which is shown in Fig. |^. 
For reference, the continuous Np is also compared with the conventional procedure 
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Figure 4: The continuous Npl^ in the analytic extension of the MS scheme as a 
function of the physical scale Q. (For reference the continuous Np is also compared 
with the conventional procedure of taking Np to be a step-function at the quark-mass 
thresholds.) 

of taking Np to be a step-function at the quark-mass thresholds. The figure shows 
clearly that there are hardly any plateaus at all for the continuous Npl^{Q) in between 
the quark masses. Thus there is really no scale below 1 TeV where Npl^{Q) can be 
approximated by a constant; for all Q below 1 TeV there is always one quark with 
mass rrii such that mf <^ or ^ is not true. We also note that if one 
would use any other scale than the BLM-scale for Npl^{Q), the result would be to 
increase the difference between the analytic Np and the standard procedure of using 
the step-function at the quark-mass thresholds. 

Figure ^ shows the relative difference between the two different solutions of the 1- 
loop renormalization group equation, i.e. (5ms((5) ~ c^msiQ)) / ciMsiQ) ■ The solutions 
have been obtained numerically starting from the world average [fl 2 8|| a-^siMz) = 
0.118. The figure shows that taking the quark masses into account in the running 
leads to effects of the order of one percent which are most especially pronounced near 
thresholds. 
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Figure 5: The solid curve shows the relative difference between the solutions to the 1- 
loop renormalization group equation using continuous Np, cXms^Q), and conventional 
discrete theta-function thresholds, Q!ms((5). The dashed (dotted) curves shows the 
same quantity but using the scale 2Q {Q/2) in Npl^^. The solutions have been 



obtained numerically starting from the world average 



0.118. 



The extension of the MS-scheme proposed here provides a coupling which is an 
analytic function of both the scale and the quark masses. The new modified coupling 
5ms(Q) inherits most of the good properties of the ay scheme, including its correct 
analytic properties as a function of the quark masses and its unambiguous scale fix- 
ing [|120|| . However, the conformal coefficients in the commensurate scale relation 
between the ay and MS schemes does not preserve one of the defining criterion of the 
potential expressed in the bare charge, namely the non-occurrence of color factors cor- 
responding to an iteration of the potential. This is probably an effect of the breaking 
of conformal invariance by the MS scheme. The breaking of conformal symmetry has 
also been observed when dimensional regularization is used as a factorization scheme 

129| , |130| | and inclusive | |131|| reactions. Thus, it does not turn 



in both exclusive H 



out to be possible to extend the modified scheme a^s beyond leading order without 
running into an intrinsic contradiction with conformal symmetry. 
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12 Application of Commensurate Scale Relations 
to the Hard QCD Pomeron 



The observation of rapidly increasing structure functions in deep inelastic scattering 
at small-Xfej and the observation of rapidly increasing diffractive processes such as 
7*p — *• pp at high energies at HERA is in agreement with the expectations of the 
BFKL 11321 QCD high-energy limit. The highest eigenvalue, u;™"^, of the leading 



order (LO) BFKL equation |p,32|| is related to the intercept of the Pomeron which in 
turn governs the high-energy asymptotics of the cross sections: a ~ = s^™""". 

The BFKL Pomeron intercept in LO turns out to be rather large: ap — 1 = uj'f^"-^ = 
12 In 2 (as/vr) ~ 0.55 for as = 0.2; hence, it is very important to know the next-to- 
leading order (NLO) corrections. 

Recently the NLO corrections to the BFKL resummation of energy logarithms 
were calculated] 133 , T^] by employing the modified minimal subtraction scheme 
(MS) [|135|] to regulate the ultraviolet divergences with arbitrary scale setting. The 



NLO corrections ||133| , |134|] to the highest eigenvalue of the BFKL equation turn out 
to be negative and even larger than the LO contribution for > 0.157. It is thus 
important to analyze the NLO BFKL resummation of energy logarithms [p.33|, p.34|| in 



physical renormalization schemes and apply the BLM-CSR method. In fact, as shown 
in a recent paper ||136|] , the reliability of QCD predictions for the intercept of the BFKL 
Pomeron at NLO when evaluated using BLM scale setting [l^ within non-Abelian 
physical schemes, such as the momentum space subtraction (MOM) scheme 111371 , |13S 



or the T-scheme based on T — > ggg decay, is significantly improved compared to the 
MS-scheme. 

The renormalization scale ambiguity problem can be resolved if one can optimize 
the choice of scales and renormalization schemes according to some sensible criteria. 



In the BLM optimal scale setting|jT6[, the renormalization scales are chosen such 
that all vacuum polarization effects from the QCD /3-function are resummed into the 
running couplings. The coefficients of the perturbative series are thus identical to 
the perturbative coefficients of the corresponding conformally invariant theory with 
(3 = 0. 

In the present case one can show that within the V-scheme (or the MS-scheme) 
the BLM procedure does not change significantly the value of the NLO coefficient 
rfz/). This can be understood since the V-scheme, as well as MS-scheme, are ad- 
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justed primarily to the case when in the LO there are dominant QED (Abehan) type 
contributions, whereas in the BFKL case there are important LO gluon-gluon (non- 
Abehan) interactions. Thus one can choose for the BFKL case the MOM-scheme 
||1371 , |138|| or the T-scheme based on T ^ ggg decay. 

Adopting ELM scale setting, the NLO BFKL eigenvalue in the MOM-scheme is 



LO 
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BLM W 1 
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(33) 
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The /3-dependent part of the VMOMi.^) defines the corresponding BLM optimal 
scale 
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(34) 



At z/ = we have Qblm'^{^) = Q^(4exp[2(l + 2J/3) - 5/3]) ~ 127. Note that 
QMEin^) contains a large factor, exp[-4Tj(^^o^^//5o] = exp[2(l + 2//3)] ^ 168, which 
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reflects a large kinematic difference between MOM- and MS- schemes [^31 

One of the striking features of this analysis is that the NLO value for the intercept 
of the BFKL Pomeron, improved by the BLM procedure, has a very weak dependence 
on the gluon virtuality Q^. This agrees with the conventional Regge-theory where one 
expects an universal intercept of the Pomeron without any Q^-dependence. The minor 
Q^-dependence obtained, on one side, provides near insensitivity of the results to the 
precise value of A, and, on the other side, leads to approximate scale and conformal 
invariance. Thus one may use conformal symmetry [|140| , |141|| for the continuation of 
the present results to the case t 7^ 0. 

The NLO corrections to the BFKL equation for the QCD Pomeron thus become 
controllable and meaningful provided one uses physical renormalization scales and 
schemes relevant to non-Abelian gauge theory. BLM optimal scale setting auto- 
matically sets the appropriate physical renormalization scale by absorbing the non- 
conformal /3-dependent coefficients. The strong renormalization scheme dependence 
of the NLO corrections to BFKL resummation then largely disappears. This is in con- 
trast to the unstable NLO results obtained in the conventional MS-scheme with arbi- 
trary choice of renormalization scale. A striking feature of the NLO BFKL Pomeron 
intercept in the BLM approach is its very weak Q^-dependence, which provides ap- 
proximate conformal invariance. 
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The new results presented here open new windows for apphcations of NLO BFKL 
resummation to high-energy phenomenology. 

Recently the L3 collaboration at LEPL3 has presented new results for the vir- 
tual photon cross section cr{^* {Q a)i* {Qb) hadrons using double tagged e"'"e^ 
e"''e~hadrons. This process provides a remarkably clean possible test of the perturba- 
tive QCD pomeron since there are no initial hadrons. [ 142|] The calculation of cr(7*7*) 
and is discussed in detail in references [|142|] . We note here some important features: 

i) for large virtualities, (7(7*7*) the longitudinal cross section a^i dominates and 
scales like l/Q^, where ~ max{Q^, Q\}. This is characteristic of the perturbative 
QCD prediction. Models based on Regge factorization (which work well in the soft- 
interaction regime dominating 77 scattering near the mass shell) would predict a 
higher power in 1/Q. 

ii) ct(7*7*) is affected by logarithmic corrections in the energy s to all orders in 
Us- As a result of the BFKL summation of these contributions, the cross section 
rises like a power in s, a oc s^. The Born approximation to this result (that is, the 
0{a1) contribution, corresponding to single gluon exchange gives a constant cross 
section, aBom oc s°. A fit to photon-photon sub-energy dependence measured by L3 
at = 91 GeV and < Q\ >=< Q\ >= 3.5 GeV^ gives ap - 1 = 0.28 ± 0.05. 
The L3 data at = 183 GeV and < Q\ >=< Q\ >= 14 GeV^ gives ap-l = 
0.40 ± 0.07 which shows a rise of the virtual photon cross section much stronger than 
single gluon or soft pomeron exchange, but it is compatible with the expectations 
from the NLO scale- and scheme-fixed BFKL predictions. It will be crucial to measure 
the Q\ and Q% scaling and polarization dependence and compare with the detailed 
predictions of PQCDM. 



13 Summary on Commensurate Scale Relations 

Commensurate scale relations have a number of attractive properties: 

1. The ratio of physical scales Qa/Qb which appears in commensurate scale rela- 
tions reflects the relative position of physical thresholds, i.e. quark anti-quark 
pair production. 

2. The functional dependence and perturbative expansion of the CSR are identical 
to those of a conformal scale-invariant theory where (3a{c(a) = and /9b(q;b) = 
0. 
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In the case of theories approaching fixed-point behavior i3a{c(a) = and (3B{aB) - 
0, the commensurate scale relation relates both the ratio of fixed point couplings 
ctA/ctB, and the ratio of scales as the fixed point is approached. 



4. Commensurate scale relations satisfy the Abelian correspondence principle 
i.e. the non- Abelian gauge theory prediction reduces to Abelian theory for 
Nc — > at fixed Cpds and fixed Np/Cp- 

5. The perturbative expansion of a commensurate scale relation has the same 
form as a conformal theory, and thus has no n\ renormalon growth arising from 



the /3-function. ||143|| It is an interesting conjecture whether the perturbative 
expansion relating observables to observable are in fact free of all n\ growth. The 
generalized Crewther relation, where the commensurate relation's perturbative 
expansion forms a geometric series to all orders, has convergent behavior. 

Virtually any perturbative QCD prediction can be written in the form of a com- 
mensurate scale relation, thus eliminating any uncertainty due to renormalization 
scheme or scale dependence. Recently it has been shown [p!44| how the commensu- 



rate scale relation between the radiative corrections to r-lepton decay and Re+e- {s) 
can be generalized and empirically tested for arbitrary r mass and nearly arbitrarily 
functional dependence of the r weak decay matrix element. 

An essential feature of the av(Q) scheme is the absence of any renormalization 
scale ambiguity, since is, by definition, the square of the physical momentum 
transfer. The ay scheme naturally takes into account quark mass thresholds, which 
is of particular phenomenological importance to QCD applications in the mass region 
close to threshold. As we have seen, commensurate scale relations provide an analytic 
extension of the conventional MS scheme in which many of the advantages of the 
ay scheme are inherited by the a-g^ scheme, but only minimal changes have to be 
made. Given the commensurate scale relation connecting a^g to ay expansions in 
5ms are effectively expansions in ay to the given order in perturbation theory at a 
corresponding commensurate scale. 

The calculation of ^y"*, the two-loop term in the Gell-Mann Low function for 
the ay scheme, with massive quarks gives for the first time a gauge invariant and 
renormalization scheme independent two-loop result for the effects of quarks masses 
in the running of the coupling. Renormalization scheme independence is achieved by 
using the pole mass definition for the "light" quarks which contribute to the scale 
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dependence of the static heavy quark potential. Thus the pole mass and the \^-scheme 
are closely connected and have to be used in conjunction to give reasonable results. 

It is interesting that the effective number of flavors in the two-loop coefficient of the 
Gell-Mann Low function in the ay scheme, l^'pl/^ becomes negative for intermediate 
values oiQ/m. This feature can be understood as anti-screening from the non-Abelian 
contributions and should be contrasted with the QED case where the effective number 
of flavors becomes larger than 1 for intermediate Q/m. For small Q/m the heavy 
quarks decouple explicitly as expected in a physical scheme, and for large Q/m the 
massless result is retained. 

The analyticity of the ay coupling can be utilized to obtain predictions for any 
perturbatively calculable observables including the finite quark mass effects associated 
with the running of the coupling. By employing the commensurate scale relation 
method, observables which have been calculated in the MS scheme can be related to 
the analytic V-scheme without any scale ambiguity. The commensurate scale relations 
provides the relation between the physical scales of two effective charges where they 
pass through a common flavor threshold. We also note the utility of the ay effective 
charge in supersymmetric and grand unified theories, particularly since the unification 
of couplings and masses would be expected to occur in terms of physical quantities 
rather than parameters defined by theoretical convention. 

As an example we have showed in Ref. |p.2CI|| how to calculate the finite quark mass 
corrections connected with the running of the coupling for the non-singlet hadronic 
width of the Z-boson compared with the standard treatment in the MS scheme. 
The analytic treatment in the V-scheme gives a simple and straightforward way of 
incorporating these effects for any observable. This should be contrasted with the MS 
scheme where higher twist corrections due to finite quark mass threshold effects have 
to be calculated separately for each observable. The V-scheme is especially suitable 
for problems where the quark masses are important such as for threshold production 
of heavy quarks and the hadronic width of the r lepton. 

It has also been shown that the NLO corrections to the BFKL equation for the 
QCD Pomeron become controllable and meaningful provided one uses physical renor- 
malization scales and schemes relevant to non-Abelian gauge theory. ELM optimal 
scale setting automatically sets the appropriate physical renormalization scale by 
absorbing the non-conformal /3-dependent coefficients. The strong renormalization 
scheme dependence of the NLO corrections to BFKL resummation then largely dis- 
appears. This is in contrast to the unstable NLO results obtained in the conven- 
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tional MS-scheme with arbitrary choice of renormahzation scale. A striking feature 
of the NLO BFKL Pomeron intercept in the BLM/CSR approach is its very weak 
Q^-dependence, which provides approximate conformal invariance. The new results 
presented here open new windows for applications of NLO BFKL resummation to 
high-energy phenomenology, particularly virtual photon-photon scattering. 
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